Magnetically Induced Thermal Rectification 
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We consider far from equilibrium heat transport in chaotic billiard chains with non-interacting charged parti- 
cles in the presence of non-uniform transverse magnetic field. If half of the chain is placed in a strong magnetic 
field, or if the strength of the magnetic field has a large gradient along the chain, heat current is shown to be 
asymmetric with respect to exchange of the temperatures of the heat baths. Thermal rectification factor can be 
arbitrarily large for sufficiently small temperature of one of the baths. 
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The problem of explaining irreversible macroscopic trans- 
port laws, such as the Fourier law of heat conduction, from 
the reversible microscopic equations of motion is one of the 
open problems in nonequilibrium statistical mechanics [1, 2]. 
The problem is still far from being settled, in particular in the 
framework of low dimensional systems, e.g. one dimensional 
particle chains, where different instances of anomalous trans- 
port can be identified due to various kinematic or dynamic 
mechanisms. 

In systems of non-interacting particles, e.g. quasi 1-d 
chaotic billiards, the validity of Fourier law has clearly been 
confirmed [3], even though such non-interacting system can- 
not exhibit local thermal equilibrium in the non-equilibrium 
steady state [4, 5]. Furthermore, interesting connection be- 
tween the Fourier law and particle diffusion has been estab- 
lished [6]. These problems are not only interesting for under- 
standing the fundamentals of statistical mechanics, but they 
may also have straightforward applications, e.g. for connect- 
ing dynamics and transport in emerging nano-technology, en- 
gineering of molecular motors, understanding and control of 
energy flow in bio-molecules, etc. 

In view of these ideas, Terraneo et al. [7] have recently pro- 
posed a mechanism for thermal rectification in an anharmon- 
ically interacting particle chain. Using an effective phonon 
approach, they have shown that a chain composed of three 
different parts may have asymmetric heat transport properties 
due to (non)matching of the effective phonon bands. This and 
related ideas have been further elaborated and improved[8, 9]. 
A further step to devise a thermal transistor has been discussed 
in [10] in terms of the negative differential thermal resistance 
observed in some anharmonic chains. 

Since these first works, other different mechanisms leading 
to thermal rectification have been described. In [11] it was 
shown that the nonlinearity of the dynamics of an asymmet- 
ric two-level system leads to an asymmetric heat flow. More 
recently in [13] it was shown that thermal rectification can 
be observed in asymmetric billiards of interacting particles. 
When the billiard is subjected to an external temperature gra- 
dient the effective interaction leads to a temperature depen- 



dence of the transmission coefficient and thus, it is possible to 
dynamically controls the transmission probability of the bil- 
liard. For this type of billiard systems rectifications as large 
as 10 3 were observed. However, the theory presented in [13] 
requires the knowledge of the microscopic transmission coef- 
ficients, which is, at best phenomenological. Another, simple 
phenomenological mechanism of thermal rectification, based 
on strong positional and temperature dependence of local ther- 
mal conductivity, has been proposed by Peyrard [12]. 

A common problem with most of these proposals is that the 
rectification factor has always been rather limited, i.e. it is 
very difficult to achieve a situation in which heat flows only in 
one direction and not in the opposite one. Moreover, the fact 
that the thermal rectification depends directly or indirectly on 
the microscopic particle-particle interaction, renders difficult 
the ability to control the power of rectification. 

In this letter we propose a novel microscopic mechanism 
for thermal rectification which works in the absence of par- 
ticle interactions. Instead, thermal rectification is controlled 
with an external non-uniform magnetic field, leading to an ar- 
bitrarily large power of rectification. In our model heat is 
carried by charged particles and the only restriction is that the 
typical mean-free-paths due to dissipative mechanisms should 
be much larger than the Larmor radius of the charged particles 
in the magnetic field. In practice this means that the tem- 
peratures of the baths should be quite low, but such regimes 
are nowadays easily accessible for example in mesoscopic 
physics, quantum dots/wires, anti-dot lattices etc. The recti- 
fication mechanism is very simple, it is based on asymmetric 
reflection of slow and fast particles off the interface between 
regions with different magnetic field intensities and predicts 
arbitrary large rectification factors for sufficiently low tem- 
peratures. 

We consider a gas of noninteracting point particles of 
mass m and electric charge e that move freely inside a two- 
dimensional billiard region. The billiard is shaped as a chain 
of equal cells as depicted in Fig. 1. Let the circular obstacles 
have radius R and centers of the discs be arranged in a hexag- 
onal array with lattice distance 4/V3, such that the billiard 



2 




FIG. 1: Geometry of the model: We consider a billiard composed 
by N symmetric unit cells, each of which is made of a rectangular 
boundary of width Sy2/V3 and length 5x = 4. Inside each cell 
there are three circulars hard scatterers of radius R disposed in a 
triangular lattice (see the text). Each cell is subjected to an external 
perpendicular uniform magnetic field of strength B, £ [0, B]. At 
the left and right boundaries the billiard is placed in contact with two 
stochastic thermal baths at different temperatures TL and Tr. The 
figure corresponds to a channel of N = 2 cells in which Bi = and 
£>2 = B. The dashed lines are drawn as a reference. 
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motion has a closed horizon for R > 1. Then we cut out a 
quasi one-dimensional billiard channel of width 5y = 2/ V3 
through the centers of two nearby rows of discs, so that one 
rectangular Sx x 5y cell of the channel, containing one half 
disk and two quarter disk obstacles, has length 5x ~ 4. Nega- 
tive curvature of the billiard boundary ensures that the motion 
in the absence of the magnetic field is completely chaotic - 
hyperbolic. 

The simplest model that we consider is composed of two 
cells (Fig. 1). The left cell contains no magnetic field, whereas 
the right cell is subjected to a perpendicular uniform magnetic 
field of strength B. In what follows we will refer as step con- 
figuration, to a channel of N cells for which the N/2 left cells 
contain no magnetic field and the N/2 right cells contain a 
magnetic field of density B. Let A denote the smallest length 
scale in the problem, namely the width of the opening between 
the neighboring cells in our model, 



A = — = - R . 



(1) 



The transmission probability of the billiard is controlled by 
the strength of the magnetic field. Consider the particles that 
cross the junction between the two cells from left to right. 
There exist a critical velocity v c : fast particles of velocity v > 
v c , always enter the right cell, and thus contribute to the left 
to right energy flow provided they reach the right end of the 
system which is coupled to a heat bath as explained below. 
Instead, slow particles of velocity v < v c , such that the gyro- 
magnetic radius g(v) = mv/(eB) is less than A/2 will be 
reflected or transmitted depending on the position at which 
they reach the interface. 

Using a statistical ensemble of trajectories, the condition 
for the critical velocity g(v c ) = A/2 can be rewritten as the 
condition giving a critical temperature 
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such that particles which are colder than T c will be reflected 
in their majority. 



FIG. 2: (Color online). Rectification index A as a function of the 
minimal temperature r m i n for a 2-cell channel in a step configura- 
tion. The maximal temperature was set to r max = 33.4275. The 
dashed line is for the no rectification value A = 1. The solid line 
corresponds to r"^ 2 - Inset: Color density plot for A as a function 
of the temperatures of both baths. Note that by construction the ma- 
trix is symmetric. For the sake of presentation the temperatures are 
rescaled by a = 1.33711. 



However, for the particles that cross the junction from right 
to left there is no condition on their velocity and they always 
enter the left cell. The above qualitative argument makes it 
clear that the transport of heat will be strongly asymmetric 
with respect to exchange of effective temperatures of particles 
on different sides of magnetic field boundary, provided the 
temperatures are strongly different, one being larger and the 
other smaller than T c . In the rest of the letter we measure 
the temperature in units of T c , and denote it as r = T/T c = 

We couple the left cell with a stochastic heat bath [3] of 
dimensionless temperature r L and the right cell with a heat 
bath of temperature t r , namely if a particle reaches the end 
boundary of the billiard that is in contact with a heat bath it 
is reflected with a velocity chosen from a distribution with 
probability densities 
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where r is the temperature of the respective heat bath in di- 
mensionless units. 

If the left bath is cold, i.e. r L <C 1, then the particles will 
most of the time remain in the left cell and there will be no 
heat current between the baths, irrespective of the temperature 
of the right bath, which we assume is larger than 1. If we ex- 
change the temperatures of the baths, then we will in general 
have some non-small heat current flowing, since cold particles 
have no problem in leaving the region with a magnetic field. 



FIG. 3: Energy density (a) and particle density (6) profiles for a 
10-cell channel in a step configuration. The bath temperatures are 



3.34277 x 10" 



33.4277. In both panels pro- 



files are shown for the positive and negative gradient for which 
tr = w (triangles) and tr = r min (circles) respectively. 



Let us measure the heat currents per particle in the steady 
state, and denote them as J + if t l < t r and J~ if the temper- 
atures are exchanged, i.e. r L > t r . Comparing the magnitude 
of these two currents we quantify the thermal rectification as 



max{|J + |,|J-|} 
min{|J+|,|J-|} 



(4) 



From our argument it is clear that rectification will be effec- 
tive if one of the temperatures is very small, say r L «l and 
the other is simply above the critical, t r > 1. It is possible to 
make a quantitative prediction on the scaling of rectification 
factor A with temperatures. When the energy current through 
the magnetic field interface is very weak, namely if t l < r R , 
the current is proportional to the transmission coefficient at 
the interface, i.e. one minus the probability of reflection. 

From our previous discussion it is clear that if r L < 1 < r R 
the particle density will be larger at the left cell (with zero 
magnetic field). This is mainly because the cold particles in 
the left cell spend long time before being able to cross the 
interface. A particle of velocity v is transmitted (not reflected) 
to the right cell if it crosses the interface at a distance from the 
upper boundary shorter than 2g(v) [14]. Therefore, invoking 
the ergodicity of the dynamics in the left cell, we can simply 
estimate the transmission coefficient t as 



2g(u) 
A 



(5) 



where q(v(t)) = ^/2mk B T min /eB and we denote by r min = 

m ™{T L ,T R }. 

However, in the reverse situation (exchanging t l and r R ) 
we have t~ ~ 1, so we can estimate the rectification 



A = t~ /t + oc 
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(6) 



In Fig. 2 we show the rectification index A as a function 
of T min for fixed value of the maximal temperature r max . We 
clearly confirm the scaling (6), indicating also that the recti- 
fication index is only very weakly depending on the maximal 
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FIG. 4: 10-cell channel in a gradient configuration: (a) Rectifica- 
tion index A as a function of the minimal temperature r m i n , with 
T m ax = 33.4275. The temperatures are given in units of T C (B), 
where B — 100 is the strength of the magnetic field at the right- 
most cell. The dashed line is for the no rectification value A = 1. 

— 1/2 

The solid line corresponds to r min . In the panels at the right the 
profiles of kinetic energy (6) and particle density (c) for the positive 
(triangles) and negative (circles) gradient are shown. 



temperature (as long as T max > 1. This can be seen in the 
inset of Fig. 2) where the rectification index A is shown as a 
function of the two temperatures of the heat baths t l and r R . 
The correctness of the scaling (6) shows that the magnetically 
induced rectification power is arbitrarily large for sufficiently 
small temperature r min . 

Even though a non-interacting system cannot reach local 
thermal equilibrium, and therefore the concept of (local) tem- 
perature is not well defined, it is instructive to compute tem- 
perature and density profiles, as long-time averages of kinetic 
energy density and particle density as a function of the hori- 
zontal coordinate along the chain. In Fig. 3 we plot the kinetic 
energy density and the particle density measured in each cell 
of a channel of 10 cells, the right half of 5 cells being in a 
uniform magnetic field. We plot the profiles for the positive 
gradient, r L < r R (triangles), and the inverted negative gradi- 
ent (circles). The positive gradient for which r L is r in , corre- 
sponds to the insulating case, i.e. to the situation in which the 
heat current is very small. In Fig. 3 we see that the insulating 
case corresponds to a large gap in the kinetic energy profile. In 
contrast, for the negative gradient the gap in the kinetic energy 
profile is much smaller and this coincides with the observation 
of a larger heat current. However, note that the energy profile 
for the negative gradient is not a monotonous function due to 
the lack of local thermal equilibrium. Moreover, the density 
profile confirms our prediction that the low current situation is 
characterized by a very small density of particles on the side 
of the magnetic field. 

We discuss now a slightly different situation in which the 
magnetic field does not change abruptly from one half of the 
system to the other, but instead changes gradually, forming a 
uniform gradient of the magnetic field. We refer to this situa- 
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FIG. 5: Dependence of the thermal rectification index A as a function 
of the size of the chain L for a chain in a .s?ep-configuration with fixed 
T min = 3.34275 x 10" 3 , r max = 3.34275 x 10 2 and fixed particle 
density, namely one particle per cell. The dashed line corresponds to 
Inset: heat current J as a function of the size of the chain L for 
the positive (triangles) and negative (circles) gradient. 



tion as the gradient configuration for which the magnetic field 
in each cell has an intensity given by B i = B(i — 1)/ (N — 1) 
for i = 1, 2, . . . , N. Note that for N = 2 the gradient and the 
step configurations coincide. 

We have performed numerical simulations for the gradient 
configuration, and we found no qualitative differences with 
respect to the step configuration. In Fig. 4 we show the de- 
pendence of the rectification index on the minimal tempera- 
ture (panel a). As for the step configuration we have found 
that for sufficiently low temperature r mln , the rectification in- 
dex again grows as A ~ 1/ ^/r min . In the panels (b) and (c) 
we show the profiles of the kinetic energy and particle density 
respectively. As expected the profiles for the gradient config- 
uration are more smooth than for the step configuration. 

From our analysis it follows that our rectification effect is 
a phenomenon which exists only in far from equilibrium sit- 
uation. In the thermodynamic limit, a vanishingly small tem- 
perature gradient is established across the system. In Fig. 5 we 
show the results of numerical simulations which indicate that, 
for fixed bath temperatures, the rectification index scales as 
A ~ 1/L. In the inset of Fig. 5 the time averaged heat cur- 
rents for the positive and negative gradients are also shown. 

Finally, it is interesting to give a quantitative estimate of 
the critical temperature in Eq. (2) in physical units. Let us 
suppose that the gas of particles inside the billiard of Fig. 1 
consists of electrons. Assuming that the dimension of the 
opening is A = lOOnm and a magnetic field of B = IT is 
applied, then the critical temperature is T c ~ 0.5K. Thus, a 



rectification power of A ~ 10 would be measurable for ther- 
mal gradient given by T min ~ 1(T 3 K and T max — 10K that 
appears accessible to nowadays experiments. 



In this letter we have presented a novel mechanism for ther- 
mal rectification. This mechanism is fairly simple. It results 
from the asymmetric behavior of the dynamics at the magnetic 
interface, that leads to a simple temperature dependence of the 
transmission coefficient. Moreover, this mechanism for ther- 
mal rectification is not based on the, macroscopically unacces- 
sible, microscopic particle interaction, but on the interaction 
with an external field, making possible an easy control of the 
power of rectification. We have shown that the thermal rectifi- 
cation power is arbitrarily large for sufficiently small temper- 
ature of one of the heat baths. 

Furthermore, the physical scales needed for optimal imple- 
mentation of our theoretical model are realizable in present 
nano-scale experiments with meso-scopic devices. We be- 
lieve that it is precisely at these scales of meso-scopic physics 
where such a thermal rectifier would find interesting appli- 
cations. This fact makes an experimental verification of the 
mechanism presented here very desirable. 
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